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Abstract 

We study a Markov random process describing a muscle molecular 
motor behavior. Every motor is either bound up with a thin filament 
or unbound. In the bound state the motor creates a force proportional 
to its displacement from the neutral position. In both states the motor 
spend an exponential time depending on the state. The thin filament 
moves at its velocity proportional to average of all displacements of 
all motors. 

We assume that the time which a motor stays at the bound state 
does not depend on its displacement. Then one can find an exact 
solution of a non-linear equation appearing in the limit of infinite 
number of the motors. 
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1 Introduction 



Recent progress in the molecular motors study brought a wave of new models 
and methods in theoretical considerations. The tools involved into considera- 
tions are spread from biochemistry and biophysics to mathematics and prob- 
ability theory. As examples, we would like to mention works [H [21 [Sj [U [6] . 
This list contains only some papers from different areas. 

We concentrate here on mathematical aspects of the problem and, more 
precisely, we describe a probabilistic model of the muscle motor which leads 
to a non-linear Markov process. The latter notion was introduced by E.P.Mc- 
Kean [7]. Comparing with the usual case, in a non-linear Markov process 
transition probabilities are related to a non-linear equation. With respect 
to the probability theory, non-linear Markov processes play main role in the 
modelling. Because of non-linearity, the model appears to have many effects 
atypical for usual Markov processes that gives possibilities either to explain 
experimentally observed properties of the motors or to predict new ones. 
Following Howard' book ^LQl we shall distinguish two classes of the molecular 
motors: processive motors and non-processive ones. The distinction between 
these classes leads to different types of the models, however differences are 
relative. Perhaps one can construct a general model yielding all features of 
both sorts of motors. We consider here a model of non-processive motors 
which concerns motors involved in the muscle activity. 

2 Model 

2.1 Informal description 

Main components of a molecular motor complex are protein molecules called 
the motors themselves which perform a motion along a molecule called fila- 
ment (or microtubule) and playing a role of railway for the motors. At last, 
there exists a cargo or back bond which play passive role of a relocatable 
ingredient. As it can be extracted from the biochemical and biophysical lit- 
erature, there are several kinds of the motor protein molecules involved in 
slightly different types of movements. 

In this work we study the motors producing muscle motions. Those motors 
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are fixed by one of tlieir end at a long protein molecule called the thick 
filament. The thick filament is fastened in a cell. Thus the motors do not 
move with respect to the cell they are located in. The second end of the 
motors may be attached or not attached to a thin filament, another long 
molecule in the cell disposed parallel to the thick filament. If both ends of 
a motor are attached then the motor can cause a tension acting on the thin 
filament if the end positions are not opposite each other. This tension is 
the cause of the thin filament motion. The binding to and unbinding out 
the thin filament of any motor are random. When binding a motor chooses 
a point on the thin filament to attach randomly. The tension the motor 
creates is defined by a displacement value between two ends of the motor. 
Therefore all bounded motors create different tensions. It is common to take 
the velocity of the thin filament moving proportional to the average of all 
bounded filament tensions. 

2.2 Formal description 

The formalization of the muscle motor construction informally described 
above can be done in the following way. We introduce a random process 
of interacted particles. Every particle represent a motor. Therefore we lo- 
cate the particles at points of Z. 

Let us start with a random process (k{t) describing the binding and unbinding 
process the particle located at A; e Z. The state space is the two-point set 
D — {0,1}, where means the unbinding and 1 means the binding particle 
state. Then the infinitesimal operator of Ck{t) is the following 2x2 matrix 



where Pr (^(t) = 1/^(0) = O) = Cbt + o{t) and Pr {Ck{t) = 0/a(0) = l) = 

Cut + o{t) . 

Then the probability Pk{t) of kh motor to be unbound at the time t satisfies 
the equation 




dPkjt) 
dt 



CbPk{t)+Cu{l-Pk{t)), 



(1) 



which gives in the steady state 



c, 



■u 



(2) 



Cb + c. 



■u 
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The state space X of kh motor consists of pairs {z, e) where 2; G M if e = 1 
and z = if e = 0. As before the parameter e indicates the bounded (if 
6 = 1) and unbounded (if e = 0) positions of the particle. We define a random 
process C,k{t) describing binding and unbinding actions and the displacement 
when binding, including the deterministic moving of a single particle. Let 
b{x) be a distribution density such that J xb{x)dx > 0. The infinitesimal 
operator of the process defining the behaviour of kth particle is 



K^)f{z,l -e)dz- f{x,£) 



(1-^)- (3) 



XX -^/(x, e)e + Cu [f{k, 1 - e) - /(O, e)] e, 
ax 

Now a constant c;, > is the rate of the particle to jump to a point of R x {1} 
from the state (0,0), that is Pr (Cfc(t) G M x {1}/Ck = (0,0)) = Cbt + o{t). A 
constant > is the rate of the particle to jump to (0, 0) from any point 
in R X {1}. That means that the rate of unbinding does not depend on the 
point X G M where the particle was attached at the unbinding moment. The 
function b{x) is the probability density of the particle to bind at the point x 
if the particle was at (0,0). Here x is the displacement of the particle with 
respect to its neutral position. Constant x is positive. It is the tension of 
the single motor molecule. It can be seen from (jS]) that when the particle 
is on R, that is £ = 1, then it is moving to the point with the velocity 
proportional to x. 

For the density pk{x,t) of the particle to be on M we obtain the following 
differential equation (see [T2] ) 

= Cbb{x)Pk(t) + X— [xpk{x, t)] - CuPkix, t), (4) 

where Pk{t) = 1 — J pk{x,t)dx. The probability Pk{t) to be unbound satisfies 

HID- 

Next we introduce the interaction between the particles. We cannot express 
the interaction in a Hamiltonian form. Instead we introduce a deterministic 
dynamic of all bound particles such that particles dynamic is highly corre- 
lated each to other. Moreover all particles are moving with the same velocity. 
To be more precise consider all particles in the interval [— A^, A] C Z and the 
configuration space il^ = X'"^'^! of all particles in [—A, A]. The space fi^r 
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is a disjoint union of the sets r^e.^v,...^^ = Jli^^'- We can consider any prob- 
ability distribution on Qjv as a collection of measures on spaces Jl£_^,...£^. 

The generator of the process involving all particles from [—N, N] is 

L[-N,N]f{ixk, ek),k = -N, N) = 

N . 

Cb / K^)f{--;{^k-i,£k-i),{z,l-^k),{xk+i,ek+i),...)dz- 

1 \T 



k=-N 



f{(xk,ek),k^-N,...,N) 



(5) 



N 



k=-N 



/(..., {xk-i,ek-i), {k, 1 - ek)ixk+i,ek+i), ...) 



-f{{xk,ek),k^-N,...,N) 



£k — 



N 



d 



vn 



Q^fi{xk,ek),k = -N,...,N)ek, 



k=-N 

where vn = ^^jW+i '^k=-N ^k^k — F. The term F in above formula means 
the velocity which an external force adds to the common velocity Vn — 
Ylk=-N ^k£k of all particles. 

Further, we use the following notations. Let (£_7v, ^at) be fixed. Then 
M° = {i : -N <i< N.Ei^O} and = (M°)^ = {i : -N <i < N,ei^ 
1}- 

Prom now we shall denote the vector ((x^, ek),k — —N, ...,N) by X. 

For a configuration X — {{xk-,Sk),k — —N,...,N) and i e define the 
configuration UiX for which the pair (xj, 1) is substituted by (0, 0); for i e M° 
and X e M define b^X as the configuration for which (0, £j — 0) is substituted 
by (x, 1). 

Let p^{X,t) be the measure density of all bound particles to be at the given 
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points at the time moment t. Then p {X,t) satisfies the equation 



N 



■ + 2^ — EiVNP [X, t) 



dt 

N 



i=-N 



Cb 



N AT „ 

J2 e^b{xi)p''{u,X, t) + c„ ^ (1 - Si) / p^'ibtX, t)dx- 

— AT A — ^ 



(6) 



i=-N 
N 



N 



C« ^ + Cb ^ (1 - Si 



-N 



i=-N 



i=-N 
N 



p^{X,t). 



Let p^ (x, t) be the probabihty density of the fc-th particle to be at the point 
X at the time moment t. 

We are interested in the behavior of the motor system for the large number 
of motors which formally corresponds to the limit N —* oo. In this limit we 
can substitute ^^at by its expectation value 



2N + 



N 

- ^ / xp^{x,t)dx - F 

i=-N 



(7) 



Plugging in v into ([6]) we get the nonlinear equation corresponding to non- 
linear Markov process describing interacting motors. 

Let us denote by 



z/x(dx) 



2N + 



1 ^ 



i=-N 



the random measure describing the motor distribution on the thin filament. 
The expectation of this random measure is called the first correlation measure 
and its density is called the first correlation function, n{x,t). It is evident 
that 



1 ^ 



(8) 



For the non-linear Markov process defined above the first correlation function 
n{x,t) satisfies the equation 



dn{x,t) ^ ^dn(pt)_ ^ ^^^^^^^^ _ ^^^^^ _ ^^^(^^^^^ 



dt 



dx 



(9) 
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where N{t) = J n{x,t)dx and v as above is 

V = —H J xn{x,t)dx + F (10) 

It follows from that N{t) and v{t) satisfy the equations 

N = Cb{l-N)-CuN (11) 
V = -k{vN + Cb(l - N)mi) + c„(F - v), 

where mi = J x6(x)dx. Evidently for t oo N and f tend to their limit 
values 



N=^^, (12) 

Cb + Cu 

_ xNrrii — F HCbrrii — F{cb + Cu) 



1 + xf Cb + c„ + xf^ 



c 



Let us study the dependence of v on Cu and Ch- It is clear that i; as 
Cu — >■ 0, and U — > F as — oo. If xmi — F > then < at = and 
hence there exists a value c° where v is negative and maximal in absolute 
value \v\ = v°. It is achieved at 



F^ t^ttla — F F 
< = \i—2+Cb^^ — (13) 



If there is no the external force F = then c° = ^c\y>c and 



(14) 



In the case F > the velocity becomes positive for large the unbound inten- 
sity c„ > 

Recall that all above values was obtained under the condition of small ex- 
ternal force F < xmi. If F > xmi then the velocity is positive at any 
C7/, . 
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Let us consider more general case when the returning force depends on x 
non-linearly, so instead of fllOp we have 

„ = /^(xM..*)dx- + F. (16) 

where p)[x) is some non-linear function. 

Now we substitute the expression f|T5l) in the equation ([9]). The expression 
for the time derivative of v is 

Lp—dx + Cb{l — N{t)) / ipb{x)dx — Cu{v — F). (16) 

Integrating by parts we obtain from the first term v J ip' {x)n{x, t)dx. Insert- 
ing the new variable w = f tp {x)n{x,t)dx we have w = f (p {x)h{x,t)dx 
and using (Q again we get an expression for w containing the integral 
/ (p {x)n{x, t)dx. 

If ip{x) is a polynomial then repeating this procedure we finally obtain a 
finite system of ordinary differential equations. 

The similar substitution is valid if ip{x) is any trigonometric polynomial or 
more generally any linear combination of quasi-polynomials (i.e. usual poly- 
nomials multiphed by sinusoidal and exponential functions [13] )• Consider 
the simplest case (p{x) = — xsin(aa;). Denoting w = J cos{ax)n{x,t)dx and 
rric = J cos{ax)b{x)dx, nis = J sm{ax)b{x)dx we have 

N =Cfc(l - iV) - CuN 

V = — xavw — xcb{l — N)'ms + Cu{F — v) (17) 

aviv — F) 
w = h Cb{l — JMjrric — CuW. 

H 

For any stationary point of this system the value is given by ( |T2l) as 
before, and v and w can be found from v = w = ^. It is interesting to find 
the number of the stationary points of (fT7|) . Also it is interesting to study 
the limit cycles for this system if they exist. 

For the case ^[x) = — xsinh(ax) the equation for A^, v, w are similar to ([Tj 
with the only difference in the sign before the term £!H(£z£1^ 
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